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Abstract 

The six-vertex model with domain wall boundary conditions (DWBC) on an N x N 
square lattice is considered. The two-point correlation function describing the probability 
of having two vertices in a given state at opposite (top and bottom) boundaries of the 
lattice is calculated. It is shown that this two-point boundary correlator is expressible 
in a very simple way in terms of the one-point boundary correlators of the model on 
N x N and (N — 1) x (N — 1) lattices. In alternating sign matrix (ASM) language 
this result implies that the doubly refined ^-enumerations of ASMs are just appropriate 
combinations of the singly refined ones. 
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1. Introduction 


The six-vertex model, introduced in [1], and solved for periodic boundary conditions 
in [2-4], see [5,6] for a review, has turned out to be of great interest also in the case where 
domain wall boundary conditions (DWBC) are imposed. These boundary conditions 
were originally introduced for the six-vertex model in the investigation of the norms of 
Bethe states [7], in the context of the quantum inverse scattering method (Q1SM) [8]. An 
important result was obtained in [9] where an exact determinant formula for the partition 
function was obtained, see also [10]. Subsequently, this result was found of fundamental 
importance in the proof of long-standing conjectures in enumerative combinatorics, due 
to the close connection of the model with alternating sign matrices (ASMs) [11-16], see 
also [17] for a review. It should be mentioned that ASM enumerations appear to be in turn 
deeply related with quantum spin chains and some loop models, via Razumov-Stroganov 
conjecture [18]; for recent results, see for instance [19-21] and references therein. 

An important information is also contained in correlation functions. However, because 
of difficulties caused by the lack of translational invariance, their computation is still an 
open problem in the six-vertex model with DWBC. Some simplifications take place when 
correlations are considered in vicinity of the boundaries [22], The simplest one-point 
boundary correlation functions were investigated in [23], where determinant representa¬ 
tions were obtained, analogous to that of papers [9,10] for the partition function. Even if 
these are almost the simplest correlations one can study for the considered model, they 
are nevertheless of interest, especially from a combinatorial point of view. 

In the present paper we pursue the investigation of the boundary correlation functions. 
Here we evaluate the two-point boundary correlation function which gives the probability 
of having particular vertex states at two specific sites in the first and last row of the 
lattice. Such a correlation function was recently discussed in [24], in the context of a 
graphical interpretation of the approach of [23]. From a combinatorial point of view this 
correlation function is closely related with the doubly refined ^-enumerations of ASMs; 
the case of doubly refined 1-enumeration of ASMs was studied previously in [25] where an 
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explicit expression for the corresponding generating function was derived. An interest in 
the doubly refined enumerations, in connection with the Razumov-Stroganov conjecture, 
was recently stressed in [26]. Here we give an explicit, general, and relatively simple 
expression for these quantities. 

More specifically, we show that the mentioned two-point boundary correlation func¬ 
tion, for generic values of the six-vertex model weights, is expressible in a very simple 
way in terms of the analogous one-point boundary correlation functions. This implies 
that by specializing the parameters of the model to the values corresponding to weighted 
enumerations of ASMs, one can directly obtain explicit formulae for the doubly refined 
weighted enumerations of ASMs from our results here. In particular, the doubly refined 
1-, 2-, and 3-enumerations of ASMs can be easily found from the corresponding singly 
refined ones. 

To derive the result we consider first the two-point boundary correlation function in 
the more general case of the inhomogeneous six-vertex model. In this case QISM can be 
applied. The homogeneous limit is performed next, and a determinant representation, 
analogous to those given in [9,10,23], is derived. Then, using standard techniques of the 
theory of orthogonal polynomials, along the lines of our recent paper [27], we obtain the 
final formula for the two-point boundary correlation function. 

2. The six-vertex model with DWBC and QISM 

In this paper we consider the six-vertex model on an N x N square lattice with the 
domain wall boundary conditions (DWBC). Recall that the six-vertex model is a model 
of arrows residing on the edges of the lattice, with the same number of incoming and 
outgoing arrows though each lattice vertex (this constraint being known as the ‘ice rule’). 
Each vertex can be in one of six possible states i — 1,..., 6, see figure [Q A Boltzmann 
weight Wi is assigned to each vertex according to its state i: the weights are usually chosen 
to obey the arrow-reversal symmetry 


wi = W 2 = a, ws = W 4 = b, w-j = wq = c. 


( 2 . 1 ) 
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Figure 1: The six states and their Boltzmann weights. 

The vertex weights are parameterized in the standard way in terms of a spectral parameter 
A and a crossing parameter rj, 

a = sin(A + 77 ), 6 = sin(A — rj), c = sin( 2 ?/). ( 2 . 2 ) 

In the case of the N x N lattice, the ice rule allows imposing domain wall boundary 
conditions to the six-vertex model. This means fixing the direction of the boundary 
arrows as follows: all arrows on the left and right boundaries are outgoing while on the 
top and bottom boundaries they are incoming. The partition function of the model, 
denoted as Z Nl is the sum over all possible arrow configurations 

Z N = a ni+n 2 b n 3 +n 4 c n5+n6 (2.3) 

where rq denotes the number of vertices of type i in a configuration and n^s satisfy 
n\ H-h n 6 = N 2 . 

To apply the quantum inverse scattering method (QISM) we shall consider the six- 
vertex model with DWBC in its inhomogeneous version, namely, when the weights of 
the vertex being at the intersection of k-th horizontal line (row) and a-th vertical line 
(column) are parameterized as 

a(A a , v k ) = sin(A a - u k + rj) b{ A a , v k ) = sin(A a - u k - ij) c(A a , v k ) = c = sin( 2 r 7 ). 

(2.4) 

The spectral parameters Ai,..., Ajv and zq,... , v N are assumed to be different within 
each set. A lattice with DWBC, and the assignment of the spectral parameters to rows 
and columns are shown in figure El Note, that we enumerate columns (labelled by Greek 
indices) from right to left and rows (labelled by Latin indices) from top to bottom. 
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Figure 2: A lattice with DWBC. 

Apparently, the partition function (and correlation functions) of the inhomogeneous 
model is a function of 2N variables, Zn = 2W({A Q }, {* 4 .}). After applying QISM the 
spectral parameters will be set equal within each set: A Q = A and u k = 0. We shall refer 
to this procedure as the homogeneous limit. 

We shall now introduce the main objects of QISM, such as L-operator and monodromy 
matrix. The L-operator of the six-vertex model is nothing but a matrix of the Boltzmann 
weights. To each vertex being intersection of ce-th column and k-th row one can associate 
the operator L a k{\ a ^k) which acts in the direct product of two vector spaces C 2 : in 
the ‘horizontal’ space HQ = C 2 (associated with the k-th row) and in the ‘vertical’ space 
V a = C 2 (associated with the ce-th column). The arrow states on the top and right edges 
of the vertex can be viewed as ‘in’ indices of the L-operator while those on the bottom 
and left edges as ‘out’ ones. Using spin up and spin down states 



as a basis in C 2 , we define L-operator by assuming further that up and right arrows 
correspond to spin up components while down and left arrows correspond to spin down 
ones. For the L-operator we have then the expression 

L ak ( A a , v k ) = sin(A a - u k + rj a z a a z k ) + sm(2r])(a~a£ + <T+a k ), ( 2 . 6 ) 
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where <j z , cr ± = (l/2)(cr x ±icr s/ ) are Pauli matrices and the subscripts in (12.611 indicate the 
space in which they act. 

The monodromy matrix is an ordered product of L-operators. We shall define it 
here as a product along a column, promoting the corresponding vertical space V a to be 
‘auxiliary’ space while the horizontal spaces IK*. will be treated as ‘quantum’ spaces; the 
space ‘K = <8> k=1 is thus the total quantum space. To define the monodromy matrix 
it is more convenient to think of L-operator as acting in V a <g) “K, moreover, writing it as 
2x2 matrix in V a with quantum operator entries acting in !K, 

(A a -u k + r/a z k ) sin(2 \rj)a k \ 

sin(2 r)) sin(A Q - v k - r) a z k ) I 

' [“] 

Here [a] indicates that this is the matrix with respect to the auxiliary space V a and 
a k ’ stand for quantum operators in TC, acting as Pauli matrices in K k and identically 
elsewhere. The monodromy matrix is defined by 

(a{ K) B{ X a )\ 

T a (\ a )= TT L «fcA a ,i/ fc = 1 ; • ( 2 . 8 ) 

11 \C(X a ) D(X a ) 

\ / [a] 

The operators H(A) = H(A; iq, ... , u N ), etc, act in TC and play a fundamental role in 
QISM. 

These operators, H(A), B( A), (7(A), and D( A) are subject to the Yang-Baxter algebra, 


Laki, Act, Vk) 


Sill 


R aa '(X, A') [T a (X) <g> T a /(A')] = [T a (A') <g> T a > (A)] R aa > (A, A'). 

generated by the six-vertex model /?-matrix, 

( A', A) 0 0 0 N 

0 <7(A 7 , A) 1 0 

0 1 ^(A 7 , A) 0 

0 0 0 /(A 7 , A)/ 

/ aa' 


Raa’{ A, A') — 


V 


where 


/(A', A) = g( X t A ) = sin ( 21 ') 


sin(A — A') 


sin (A — A') 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 
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Figure 3: Graphical interpretation of the operators AL(A), 5(A), (7(A), and 5(A). 

The relation (USD is also known as RTT relation and it is a consequence of the similar RLL 
relation. Among the sixteen relations contained in m, in particular, are the following 


5(A) 5(A') = 5(A') 5(A), (2.12) 

A(A) 5(A) = /(A, A') 5(A) A(X) + g( A, A) 5(A) A(A'), (2.13) 

5(A) 5(A') = /(A', A) 5(A') 5(A) + g( A, A') 5(A) 5(A'), (2.14) 

C'(A) 5(A) = 5(A) 5(A) + g{ A, A) [A(A) 5(A) - A(A) 5(A)], (2.15) 

which will be used below. 


The operators A(A), 5(A), (7(A), and 5(A) admit simple graphical interpretation 
as columns of the lattice, with top and bottom arrows fixed, see figure 01 Taking into 
account this nice graphical interpretation and exploiting further the correspondence be¬ 
tween arrows and spin states to express arrows on the left and right boundaries of the 
lattice in terms of ‘all spins down’ and ‘all spins up’ states, 

N N 

<-H= ® fe(l I, It) = <8> I T )fc, (2.16) 

k =1 k =1 

where If)*, and | J, )fe are the spin up and spin down vectors of the space “Kk, it can be 
easily seen that the partition function, Z N = Z N ({ A Q };{z/ fc }), of the six-vertex model 
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with DWBC is given by 




(2.17) 


Note, that because of relation (12.121) the actual order of operators B{ A) is not important 
in this formula. 

The explicit expression for the partition function, found and proven in [9], see also [10], 
reads 

rio=i nJLi a ( A «> vk)K\ a , Vk) 


Z N = 


ni<a</3<jV A o) ril<j<fe<JV ^( U ji V k) 


det n T 


(2.18) 


where 


d(A, A') := sin(A - A') (2.19) 

and the functions a(A, v) and b( A, v) are defined in (12.21) . The N x N matrix T is given 
by 


T ak = t(X a , I/fc), t( A, v) = — 


sin (2 rf) 


■ ( \ M \ -V ( 2 ' 2Q ) 

sm(A — v + rj) sin (A — zv — 77) 

In [9,10] the determinant formula (12.1 81) was proven on the basis of certain recursion 
formulae, established in [7]. Another proof, solely on the basis of Yang-Baxter algebra, 
was given in [23]; in the next section we sketch the derivation. 

I 11 the homogenous limit, when X a = A and Vk = 0, expression dSHHD turns into the 
following one [10] 

kin IX — rd sin ( \ -A- rdl ^ 2 

( 2 . 21 ) 


^ [sin(A — rf) sin(A + rj)] N 

- - T-rAf-1, , no - aet N <P 


n ;;=>!) 5 


where 


®ak = dt +k 2 <^(A, 77 ), </?(A, rj) — 


sin (2r?) 


( 2 . 22 ) 


sin(A — rf) sin(A + rf) 

The procedure of obtaining (ESP from (ITTH1) was explained in detail in [10]. In our 
treatment of the correlation functions we shall proceed in the same way, first obtaining 
an expression for the inhomogeneous model, and next taking the homogenous limit. 


3. One-point boundary correlation functions 

Here we recall the main results of paper [23] on one-point boundary correlation functions. 
In the next section we shall explain how this approach can be used to compute two-point 









boundary correlation functions. 

In paper [23] two closely related kinds of one-point boundary correlation functions 
were considered. The first correlation function, denoted as H N \ is the probability of 
finding the r-th vertex (counted from the right) on the first row in the state i — 5. 
Formally, this correlation function can be defined by 

HV = ZJ } 1 ^\B(X n ) ... B( A r+1 ) qi B(X r ) pi B(X r -i) • ■ ■ B( A,) | ft ) (3.1) 

where qj. and pk denote projection operators on the spin up and spin down states, at k- th 
‘site’, respectively, 

% = \( 1 -Vk)i Pk = \( 1 + °k)- ( 3 - 2 ) 

(V) 

The second correlation, denoted as G N , is the boundary ‘polarization’, i.e. the probability 
of finding an arrow pointing left on the horizontal edge of the first row between r-th and 
(r + l)-th columns. One can define this correlation function by 

GV = Z- N \^\B{X N )... B{ A r+1 ) qi B(X r ) • ■ • B(Xi )\t >■ (3.3) 

Due to DWBC, the two correlation function and G$ are related by 

r 

z _ 'f( r ) \ ^ rr( a ). tr( r ) (o /|\ 

(j-jY — 2_^ n N ) n N ~ 

a=l 

These relations can be easily found by exploring the graphical interpretation of these 
functions, see figure 0 

To compute the one-point boundary correlation functions one can introduce the op¬ 
erators A(A), B( A), C(A) and D{ A) as entries of the monodromy matrix on N — 1 sites, 


V 

Tq(Aq,) . L a k (A q , z//j) 

( ^4(A a ) 

B( X a )\ 

(3.5) 

k =2 

\C(A q ) 

D(X a )J ' 


Correspondingly, one can define 




~ N 

~ N 



(tl = ® *<TI, 

It) = ® | 

T )k ■ 

(3.6) 

k =2 

k =2 



One has 

N 


N 


-4(A)! t) = I[a(A^)|5), 

5(A)|t) = 

6(A, iz fc )| ff). 

(3.7) 


k =2 k =2 
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1 


r ■<- 


1 



Figure 4: Boundary one-point correlation functions: (a) function H^p ; (b) function . 

In what follows, to simplify writing, we shall omit tildes over operators and vectors. 

The role of these new operators in application to the one-point boundary correlators 
is quite clear from their graphical interpretation shown on figure |3] and the analogous 
graphical representation for the correlation functions shown on figure 0 Specifically, in 
the case of HP all vertices of the first row are fixed, thus leading to the following formula 
in terms of the operators on N — 1 sites 

N r —1 

hP = Zp c a(X a , vi) b(X a , ui) 

a=r +1 a=l 

x (^ \B{X n ) • • • B(A r+1 ) A(X r ) B(A r _0 • • -B( A,) | ft ). (3.8) 

(V) 

Similarly, the correlation function G N can be written as a sum of r such terms 

r N a—1 

a=l 0=ct+l /3=1 

x (^ \B(X n ) ■ ■ ■ B(X a+1 ) A(X a ) B( A q _0 • • • B(Ai)| ft ). (3.9) 

On the hrst stage of computation the correlation function are expressed in terms of 
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the partition functions on (N — 1) x (IV — 1) sublattices. Using the commutation relation 
m we have 

r-l r N (\ \ \ T r 

mk ) n B (v)i fr >= tt i n /< A «. a h) n b m \ n >• p-w) 

a=l a=lk=2 ai r ’ 0=1 0=1 

P+ol / 3 ^ot 

(r) 

The last relation allows one to obtain, in the case of H N , the expression 


N r —1 

hP = Z N l c JJ a(A Q ,z/i) Jj6(A a ,i/i) 

a=r +1 a=l 

x X^ n/(V A 7 )Z iY _i({A 5 }^i )< 5 7 4 Q! ; {^}f =2 ) (3.11) 

p = lk =2 r ' 7=1 

l^p 

The corresponding expression for Gp can be immediately obtained using the relationship 
with the function H ( p , see (Id.411 . However, as explained in [23], in this case one gets the 
result in terms of some double sum which is not actually an analogue of (EH). It was 
pointed out that the analogue of (ETTT1) can be found if one uses instead the relation 


r r a—1 

x n a(Xg, u i) K^P ) u i) B(X r ) ■ ■ ■ B(X a+1 )A(X a )B(X a _ 1 ) ■ ■ ■ .B(Ai)| ff) 

ql= 1 (3=a+ 1 (3=1 

r N r r r 

=xn Aq ., Vk ) 11 b ( Xp , ) n/(A..v)n^)i« < 3 - 12 ) 

a=l k =2 /3=1 /3=1 /3=1 

P^a P^ a P 7 ^ q 

which can be obtained by means of by taking into account that LHS here is 

symmetric with respect to permutations of Ai,..., A r by construction (for more details 
on derivation of such relations see, e.g., [8], sections VII.2 and XXII.2). In this way one 
obtains 


N r 

G { P = Zp J] a(X a , vi) H&(A a ,z/i) 

a=r+1 ot=1 

r N r 

x XI J7 \——r n a ( A V ^ II/(V>-^ 7)^-1 <5^/?; {^ 1 ^= 2 )- (3.13) 

/3=1 /3 ’ 12 fc=2 7=1 

77^/3 

As explained in paper [23] this representation is important since it allows one to prove 
the determinant formula (12.1811 . Indeed, since by definition G ^ = 1, expression (13.131) 
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turns into some recurrence relation connecting the partition functions Z^ and Z^-i, 
note, that here all the parameters in the sets {A„}^ =1 and {v k }f? =1 are assumed to be 
completely arbitrary (cf. [9,10]). As explained in [23], it can be readily shown, in virtue 
of the Kramer rule and of some particular identity, that the solution of this recurrence 
relation, with the initial condition Z\ = c, is given by determinant formula (irm 

Thus, the expression for the partition function being proved within the considered 
framework, it can be used to obtain similar representations for one-point functions. Sub¬ 
stituting the expression for Z N _i in (13.1111 and (13.1311 gives rise to the determinant formu¬ 
lae for one-point functions, see again [23] for details. We end up this section by quoting 
the results. 

The function is given by 



(3.14) 


where the matrix V differs from the matrix T, equation (1 TM . just by the elements of 
the first column, 



(3.15) 


Here the function u r (A) is given by 



(3.16) 


where 


e(A, A') = sin(A — A' + 2 rj) 


(3.17) 


and the function d( A, A') is defined in (12.191) . 

(r) 

The function G N is given by 



(3.18) 


where the matrix S also differs from T just by the elements of the first column, 











Here the function s r (A) is given by 


n«„. +1 <i(A„,A)n«, 1 e(A„,A) 


s r (A) 


nii Vk) 


(3.20) 


(r) 

In the homogeneous limit the following formulae are valid. The function H N is given 


by 


H {r) ~ 

n N ~ 


(N — 1)! sin( 2 ? 7 ) 


detN& 


(3.21) 


[sin(A + 77 )] r [sin(A — 77 )]^ r+1 detjv$ 
where the matrix & differs from the matrix <P, equation (E 22 ). just by the elements of 
the last column 

a _ 1 (sin£) Ar_r [sin(£ — 2//)] r_1 


^ a,k = $a,ki k — 1 , . . . , N -* 1 ; ^a,N = d' t 


(r) 

Similarly, the function G N is given by 


M r ) _ 
U !\ ! ~ 


(JV-l)! 


[sin(e + A — 


det/v<9 


. (3.22) 


£=0 


[sin(A + r])] T [sin(A — 77 )] 


N—r 


det)v$ 


(3.23) 


where 


&a,k = ^ a,k, k = 1,..., N — 1; 0 a ,n = —91 


1 (sin£) iV r [sin(£ — 2r])) r 


[sin(e + A - 77 )] 


N 


. (3.24) 


£=0 


4. Two-point boundary correlations 

The one-point boundary correlation functions just considered can be directly generalized 
to the case of two-point ones [24], Here we shall consider in detail derivation of the 
function which gives the probability of hireling vertices of type i — 5 on the 

opposite, top and bottom, boundaries. More precisely, we define i/^ 1,r2 ' ) as the probability 
to find vertices of type i = 5 both at ri-th position (counted from the right) of the first 
row and at r 2 -th position of the last row, see figure [21 Formally, if rq < r 2 then one can 
define this correlation function by 

H^’ r2) = I B(X n ) • • • B(X r2+1 ) q N B{ A r2 ) p N B(X r2 _ l ) • • • 

x ■■■B(X ri+ 1 )q 1 B(X ri )p 1 B(X ri _ 1 )---B(X 1 )\it). (4.1) 
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Figure 5: The two-point correlation function H^ 1,r2 \ 

For r\ = r 2 or r\ > r 2 one has analogous expressions. 

Quite similarly, one can also consider the function G^ 1,r2 \ giving the probability of 
finding arrows pointing left on the first and the last rows. The functions and 

are related to each other just as in the one-point case; here we focus on the 
function H^ 1 ’’’ 2 ^ since it is the most interesting for combinatorial applications. The result 
for G^ 1 ’ r2 ' > can be found using its obvious relation with H^ 1,r2 \ 

To compute the function H ^ 1,r2 ' > we introduce the operators AL(A), B{ A), C(A) and 
D( A) as entries of the monodromy matrix on N — 2 sites, 


JV-1 

Ta(\a) • II L a l^k) 

_ /A(A q ) 

5(A Q )\ 

(4.2) 

k=2 

\C(A„) 

5(A Q )y' 


Correspondingly, we define 




~ JV —1 

~ JV —1 



(-IM = ® fc(TI, 

|t} = ® 

IT)*- 

(4.3) 

k=2 

k=2 



We have 

AT—1 


7V_1 


-4(A)| t) = a(A,i/ fc )| t), 

D( A)|t) = 


(4.4) 


k =2 k =2 
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In what follows, to simplify writing, we shall again omit tildes over operators and vectors. 

On the first stage of computation we express the correlation function H i ^ 1,r2 " > in terms 
of partition functions on (TV — 2) x (TV — 2) lattices. This can be done by directly 
following the procedure of the previous section. For an alternative way of performing this 
preliminary step, see also [24]. 

Let us assume that ry 7 ^ r 2 . If ry < r 2 then, using the same arguments as for one- 
point correlations (e.g., by invoking the graphical interpretation of the operators shown 
on figure 02 ), for the two-point correlator we write 

N 7*1 — 1 N V2 — 1 

Hp ,r2> = Z^ 1 c 2 JJ a(A Q ,z/i) JJ b(X a ,u 1 ) JJ 6(A a , u N ) JJ a(A a , u N ) 

a=n+1 a= 1 a=r2+l a=l 

x {4 |B(A„) ■ ■ ■ B (\„ +1 ) D(\„) B( A„_j).. ■ B( A ri+1 ) Al(A„) B(A„_i) ■ ■ ■ B( A,) ft). 

(4.5) 


Acting hrst with the operator A(A ri ) on the right, using 

r —1 r N—l /, , \ r 

9\^on ^r 


A(A r ) n b(a„)i d >=x n a < A «^ , v t n /<v, v> n b (v>ih->. (4.6) 

a=l a=l k= 2 a, r) o =1 fl=l 


A=1 

0 +ol 


0=1 

0 ^a 


and next with the operator TA(A r2 ), using 

r —1 r N—l 


dm n B (v)i =x n 6 ( a « vt) 3 f(xM n /(a * a ») n b <v>i * >. < 47 > 


a=l 


a=l k =2 


/9=1 


A=i 


we straightforwardly obtain 


iV 


Hp ,n) = Z N X c 2 JJ a(A Q ,i/i) JJ 6(A a ,i/!) JJ 6(A a , zyv) JJ a(A a , zyv) 


ri —1 


AT 


7*2 — 1 


a=r i+l 


ri r2 AT—1 


Q^^ort ^ri) 


Q!—1 

r 1 


cx.=r 2 +1 
AT-1 


CH—1 


x EEII a ( A “’ z/fc ) / ( A A ) 11 

OL= 1 / 3=1 fc =2 7 v ri ' 7=1 


n/< a » a j n h. v. n /(v,^ 

7=1 fc=2 J\ rn 0) 7=1 


r 2 


0 ^a 


1+OL 


7 = 
l¥=oi ,0 


x ^({A,}^,^;!^}^ 1 )- ( 4 - 8 ) 


If ri > r 2 then one should hrst act with the operator TA(A r2 ), and next with the operator 
x4(A ri ). It can be easily verified that one obtains exactly the same expression, i.e., formula 
(Oil is valid for ry 7 ^ r 2 . 
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Consider now the case r\ = r 2 = r. In this case instead of ESI) we have 

TV r —1 

H^' r) = Zjf 1 c 2 JJ a(X a , Vi) b(X a , v N ) JJ 6(A a , v x ) a( A a , v N ) 

a=r -\-1 ol— 1 


x (^ \B(X n ) • ■ • B( A r+1 ) C{K) B(Xr-i) ■ • • B(Xi )| ft). (4.9) 


Using commutation relation (IXT51) the following formula can be derived (see [8], section 
VII.2.2) 


r— 1 r TV—1 TV—1 

C(A r )[]fi(A a )|t) = eii fl(Aai Vk) b(Xfi, Vk) 

a=l a,/3=l fc=2 fc=2 

x 7tHHt7(Hi4/(v,a^) n /(a„, v)/(v,Arf n (4.io> 

J(A q , A r )/(,A r , ^ =1 7=1 

l¥=ot,l3 1^01,13 

Thus, in this case we have 


IV 


r—1 


= Zjy 1 c 2 a(A a ,i/i)6(A Q ,i/iv) X a (-V, ^jv) 


a=r+l 
r V—1 


a=l 


V—1 


X 


En a(A a , Vk) ^(A^, Vk) 


a, /3=1 fc=2 

a#/3 


k=2 


fi , (A q, A r )g(A r , A^) 
/(Aq,, A r )/(A r , A/?) 


/(Aa, A/3) 


r 

X n /(A«,X > )/(A 7 ,A fl )Z i ,_ 2 ({A { }t I ^;{^}Ji 1 ). (4.11) 


7=1 

7^a,/3 


This formula shows that expression ESI) is, in fact, valid for all values of r i and r 2 , 


1 < n,r 2 < V with no further restriction. 

Substituting now the determinant representation 112 .1 81) for the partition function Zn 
and that for the reduced partition function Ejy_ 2 into ESJh and cancelling whatever 
possible, we find 


N -1 


H 


(r 1 , 7 - 2 ) 
N 


c 2 d{v 1 , 1 /jv) n d(vi,v k )d(v k ,v N ) 
k =2 


rp TV r2 TV 

n a(A Q ,z/i) n K^a,vi) n K^^n) n a(A Q ,z/Ar) detjvT 

ct=l a=r 1 a=l a=r 2 


X 


illV+a+/3„ 3V ri (X Q )w r2 (Xp) fA -,n\ 

e «/3-- detjv-22( a ,/3 ; i,Ar). (4.12) 


«=1 /3=1 
/¥« 
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Here Tt a ,/ 3 ;i,N) denotes the (N — 2) x (N — 2) matrix obtained from T by eliminating the 
a-th and /3-th rows and the first and the last columns (i.e., its determinant is just the 
corresponding minor of T). The symbol e Q( a stands for the sign function 


1 if a > (3 
e ap — 0 if a = (3 ■ 

-1 if a < j3 

The functions w r (X) and w(X) are defined as (compare with (18.161) ' 

m nL r+ i4A«.A)nir=ho«,i) 

w r (\) = - 


(4.13) 


w r (X) = 


Yl^KKvk) 

nOlr+i lO. A a ) I e(A, A„) 


nh ! a ( A - " i . i 


(4.14) 


(4.15) 


Note, that both w r (X a ) and w r (X a ) vanish if r < a; thus the double sum in (gZEZD can be 
formally extended to N. 

Consider now how the homogeneous limit can be performed in the just obtained 
expression for the two-point correlation function. Here it should be mentioned that, con¬ 
trary to the formulae for the partition function and one-point correlators, the two-point 
boundary correlation function cannot be expressed as a determinant in the inhomoge¬ 
neous model. Nevertheless, as we shall show now this is just a minor inconvenience which 
can be very efficiently solved without much efforts. 

First of all we mention that the problem we are actually facing is that of finding the 
homogeneous limit for the following double sum 


N 


w = (-1 ) N+a+p e a/3 F(X a ,X^) det N _ 2 T {aAhN) 

a,0=1 

where F( A, A') is some function. Consider the substitution 

Aq, — A T z a , ol — 1,..., N. 


(4.16) 


(4.17) 


In what follows we shall send all z a (and Uk) to zero but before performing the limit we 
note that the double sum above can be rewritten, using 


F(X + zi, X + z 2 ) = exp (z a d £l ) exp(zpd £2 )F(X + e 1} X + e 2 ) 


£ 1=0 ’ 
£ 2=0 


(4.18) 
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in a quite formal way as follows 



exp(£i<9 £1 ) 

t( A + zi, 12 2 ) 

■ t(\ + Z!, 12 N _!) 

exp(x,<9 £2 ) 

w = 

exp(z 2 <9 £1 ) 

t(X + z 2 , v 2 ) 

t(A + Z 2 , 12N — 1 ) 

exp(^ 2 a £2 ) 


exp (z N d £l ) 

t( A + z N , 12 2 ) . 

. f(A + zn, i 2 n-i) 

exp(^jv<9 £2 ) 


i^(A + £i, A + £ 2 ) 


£ 1=0 

£ 2=0 


(4.19) 

so that the homogenous limit of expression dHI2D can be performed exactly along the 
lines given in [10]. 

Namely, we shall put z/ fc ’s and z& s to zero in the order zq, ..., i2 N , zi ,..., Zn each time 
keeping the leading order. The prefactor in (gZEZD becomes 


sin 2 (2//) 


[sin(A + r 7 )] iV+r ' 1_r2+1 [sin(A — rj )] N+r2 - ri+1 det N @ 


x 


(N — 1)! (N — 2)! 


^2 ' 


z 2 z 3 
z 3 ~4 

2! 3! 


(N- 1)! 


Z / 3 • 


JV-2 
V N-1 
(TV-3)! 


(4.20) 


while the double sum goes into 


£S . £4 ... 


v N-l 

Z N 


2! 3! (N — 1)! 


Va 


"3 ■ wr ■ • • 


,.N -2 
AV-1 


2! (N — 3)! 


x 


det (<P, 


a,k 


d: 


\Ol—l 
£2 


a: 


a— 1 
£1 


l<a<N,l<k<N—2 


7,<>l,r 2 ) 

n N 


(Si, £2) 


(4.21) 


where 

h^ 1 ' £ 2 ) — 


(sine^ ri [sin(e 1 — 2 r/)] ri 1 (sine 2 ) iV r2 [sin(e 2 + 2-//)] 


£l=£2=0 


r 2 —1 


sin(e 2 - £1 + 2r7)[sin(ei + A - T7)] Ar_2 [sin(£ 2 + A + r ])} N ~ 2 ' 
Thus, in the homogeneous limit for the two-point correlation function we obtain 


(4.22) 


H 


( ri,r 2 ) 


(N - 1)1 (N - 2)1 sin 2 (2 rj) 


N 


N+r i—r 2 +l 


[sin(A + rj)] 

det ( <P a i. 


[sin(A — r /)] 


A r +r 2 —ri+l 


det^ 


x 


<9: 


\a —1 


o: 


, 0-1 


h 


{ri,r 2 ) 

N 


( £ 1) £ 2 ) 


(4.23) 


£ 1 =£ 2=0 


l<a<N,l<k<N—2 

This determinant representation is analogous to those of the previous section for one-point 
correlation functions. 
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5. Orthogonal polynomials representation 


In this section the results for one- and two-point boundary correlation functions will be 
analyzed by making use of the orthogonal polynomials theory, along the lines proposed 
in paper [27]. Here we show that the two-point boundary correlation function, studied in 
the previous section, is expressible in terms of one-point ones. 

Let {.P„(:c)}^L 0 be a set of polynomials, with non-vanishing leading coefficient 

P n {x) = K n x n + ..., kappa n ^ 0, (5.1) 


and orthogonal on the real axis with respect to some weight p(x), 

/ OO 

P ni (x) P+x) /jj(x) da: = h ni S nin2 . 

-OO 

Let c n denote n-tli moment of the weight p(x), i.e. 

/ OO 

x n /j,(x)dx , n — 0,1,... 

-OO 

and let us consider the (n + 1) x (n + 1) determinant 


(5.2) 


(5.3) 



Co 

Cl • 

• c n 

= 

Cl 

C 2 

• c n -\-i 


Cn 

Cn+1 

C 2 n 


(5.4) 


Using the orthogonality condition (El and well-known properties of determinants, one 
can easily find the following formula 


hi 


A, = IB- 


k =0 K 


(5.5) 


This formula can be used for computation of determinants, provided the orthogonal 
polynomials {P n (x)}™ =0 are known. On the other hand, the polynomials {P n (x)}^L 0 can 
in turn be expressed as determinants. For later use let us introduce the notation 


^n\x i,...,X k ) = 


Co 

Cl . 

• C n —\ z 

l 

1 . 

. 1 

Cl 

c 2 • 

• C n — ^_|_i 

X\ 

X 2 • 

• X k 

Cn 

Cn+1 • 

• C2n—k 

rv^n 

X 1 

r+n 

ry>n 

■ x k 


(5.6) 
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(5.7) 


so that An = A n . For the polynomials one can find that 


Pn{x) = -p- 

^n -1 


A£\x). 


For a proof, see, e.g., book [29]. 

The relation El can be read off inversely thus giving an expression for the deter¬ 
minant An\x) in terms of the polynomials P n {x). Taking into account that (see El) 


A r . 


A n -1 


we can write 


A (1) 


n \X) K n . . 
= T P n{x) 
iln 


A 


(5.8) 


(5.9) 


( e Y\ 

Consider now the case of An [xi,x 2 ). It is clear that the term of the highest powers 
on both x\ and x 2 is just A n - 2 (x "x] 1-1 — extending further the methods of [29] 

we obtain 


A^( Xl ,x 2 ) = 


A n -2 

^n^n —1 


[f’n-l(ll)7’,.fe) - Pn(x,)P n _!(SC 2 )] 


(5.10) 


Again using El , we write 

An\x\,X 2 ) _ K n K n -i 

— h h 
1 i n li n—l 


A r 


^n^n —1 

hnhn —1 


[.P n -l(xi)P n (x 2 ) - Pnix^Pn-!^)] 

Pn—l{Xl) Pni^X l) 

Pn- l(x 2 ) PnM 


(5.11) 


This formula can be easily extended to the general case of A^ n \x i,...,x n ); in what 
follows we shall make use only of formulae El and EH- 

Consider now how we can use all these formulae in application to the boundary cor¬ 
relation functions. First we note, following paper [27], that the determinant entering the 
expression for the homogenous model partition function can be related with orthogonal 
polynomials using the integral representation 


sin (27/) 


c(A—77/2) sinh(Tyx) ^ 

sinh(7rx/2) 


(5.12) 


sin(A — rj) sin(A + rj) J_ 

This formula is valid if 0 < rj < 7t/ 2 and ?/ < A < tt — rj] these values of A and rj 
correspond to the so-called disordered regime of the six-vertex model (for similar formulae 
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valid for other regimes, see [28]). This regime is the most interesting especially for 
combinatorial applications of the six-vertex model with DWBC (in these applications 
one further specializes to A = 7 t/ 2). It can be easily seen that our results below, however, 
do not depend on the particular choice of the regime, and can be extended to other 
regimes simply using the proper analytical continuation in the parameters A and 77 . 

Formula (EH implies that we have to deal with the set of polynomials which are 
orthogonal with respect to the following weight function 


H(x) = n{x\\,7}) 


MX-n/ 2 ) sinhfos) 
sinh(7nr/2) 


(5.13) 


The corresponding polynomials P n {x) = P n (x; \rj) also depend on A and rj which are to 
be considered as parameters. For later use let us mention the following useful property 
of these polynomials 

P n (x; A, rf) = (-l) n P n {~x\ 7 r — A, 77 ). (5.14) 


This property can be easily established in virtue of formula m- It is to mentioned 
also that both the leading coefficient K n = K n (X,rj) and the normalization constant h n = 
h n ( A, r/) are invariant under the substitution A —> n — A. 

The transformation A —> n — A is related to the so-called crossing symmetry of the 
six-vertex model which has useful consequences for the one-point boundary correlation 

(r) 

function H N . Recall that the the crossing symmetry is the symmetry of the vertex 
weights under reflection with respect to the vertical axis, and simultaneous interchange 
of the functions a and b, which is equivalent to setting A —> tt — A. As we shall explain 
now, these simple properties related to the crossing symmetry allow one to derive easily 
two equivalent representations for the one-point boundary correlation function H N . Us¬ 
ing these formulae we shall then show that the two-point boundary correlation function 

1,r2 ' ) is expressible in terms of one-point boundary correlators. 

Indeed, since the lattice with DWBC is invariant under the reflection with respect to 
the vertical axis, the crossing symmetry thus imply that the following relation holds 


hP( a, if) 


(jV-r+l) 

n N 


(7T 


A, 77 ). 


(5.15) 
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Consider expression (Id. 21 1 ) for the one-point function H^\ Due to (15.91) we can rewrite 
it as 


fffdvf) = 


(N — 1)! sin(2?/) 


«jv-i(A,T7) 


[sin(A + r])] r [sin(A — r])] N r+1 ^n- i(A, 77 ) 


x P/v_i(< 9 e ; A,77) 


sine 


\N—r 


sin(e — 2 77 )] 


r—1 


[sin(e + A - rjj] 


N -1 


(5.16) 


e=0 


Taking into account ( 1011 ) and the properties of the leading coefficient n n {\ , 77 ) and the 
normalization constant h n {X,if mentioned above, it can be easily seen that from (15 .1 51) 
and (15.161) the following expression is valid as well 


1 $ ) M = 


(N — 1)! sin (2if) 


K N -i(\,T j) 


[sin(A + r])] r [sin(A — r])] N r+1 ^iv-i(A, 77 ) 


x P N _ i(9 e ; A,?/) 


sm e 


\T — 1 


sin(e + 277 )] 


N-r 


[sin(e + A + r))] 1 *- 1 


(5.17) 


£=0 


Note that this expression means simply that the limit e —> 0 in (15.161) can be changed 
into £ —> 277 without altering the result. 

Thus, these two equivalent representations, (15.161) and (1071) . can be used in the study 
of the two-point correlation function H^ 1,r2 ^ given by expression (|4.2dl) . which certainly 
involves similar structures. 


Before turning to this analysis, let us put the above formulae for the one-point function 
in a more compact and convenient notations. In what follows we shall often omit the 
dependence on A and 17 where possible. 

We define the functions 


sin(A + if sine sin(A — 77 ) sin(e — 2 ?/) 

sin(A — if sin(e — 2 rf) ’ sin( 2 ? 7 ) sin(e + A — if)' 

which are related to each other as 

Also we define 


(5.18) 


(5.19) 


sin(A — rf) sine . ^ sin(A + rf) sin(e + 2 rf) 

sin(A + rf) sin(e + 2 rf) ’ sin( 2 ? 7 ) sin(e + A + 17 ) ’ 
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which are in turn related to each other as 


q(e) = 


1 -&(e)' 

Note, that the functions with tildes are introduced such that 


(5.21) 


w(e;A,r;) A, if), g(e; A, >j) = —g(-e; ir - A, if) (5.22) 

in accordance with the crossing symmetry considerations made above. Additionally, let 
us denote 

Kn-i(x) = (N- 1)! ^ ^ Pjv-i(z) (5.23) 

tlN-l 

where ip = p(X,7]) is exactly the function defining entries of the matrix A, see (12.221b In 
these notations formulae (15.1611 and (EUD for the correlation function read 

H { ; ] = (5.24) 

and 


£=0 


HP = K N . 1 {d e )[u(e)] r L [g(e)] 


r —1 r 


1 AT—1 


£—0 


(5.25) 


respectively. 

Consider now our main object, the two-point correlation function which is 

given by formula (jPHp . Obviously, function @22) contains all the structures introduced 
above apart from the factor sin(£ 2 — £i + 2r?) standing in the denominator there. However, 
using the identity 


sin( 2 ? 7 ) sin(6 2 — £i + 2 r/) = sin Ei sine 2 — sin(£: 1 — 2 r]) sin(e 2 + 2 rf) 
it can be easily seen that 

8111 ( 6 ! + A — 7?) sin(e 2 + A + 7?) 1 1 


(5.26) 


(5.27) 


sill(£ 2 - El + 2 l /) tpg(e i)0(e 2 ) - 1 

Thus, taking into account formula we can write the two-point correlation function 

in the form 


H 


(ri,r 2 ) 


N 


[A^-i(a £l )A^_ 2 (a £2 ) - Kn-2 (<9 £i ) K N -i(d £2 )] 


N—2 


X 


o;(6i)0 (6 2 ) 


. (5.28) 


£ 1 = 0 ,£ 2=0 
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Now taking into account that u(£),u(e) —> 0 as £ —> 0 we can expand the denominator in 
(EH) in power series and it can be easily seen that only the first few terms (actually not 
more than N) of this expansion will contribute. As a result, in virtue of relations 115 .1 91) 
and flEHJ) , we arrive to the following expression in terms of the one-point functions 

N 

TT(ri,r 2 ) _ ( Tr(ri-j+l) Tj-(N-r 2 +j) _ Tr{ri~j) Tr(N-r 2 +j) 

11 N ~ \ n N n N~\ n N n N -1 

j =1 

- H^_- j) H^~ r2+j+1) + (5.29) 

where it is assumed that ifr<0orr>A^ + l then Hjp = 0 by definition. The formula 
is our main result here. 

As a comment to this result let us rewrite it in terms of the generating functions. Let 
us introduce the generating function for the two-point correlation function 

N 

H n (u,v):= J2 (5.30) 

r=l,s=l 


and dehne the generating function for the one-point correlation function 


N 


h n (u) --=J2 h 


(N—r+l) r —l 

N u 


(5.31) 


r=1 


Then 1)5.291) implies that 


H n (u,v) = 


(u - 1 )H n (u) ■ vH n _ i(u) - uH n _i(u) • (v - 1)H n (v) 

u — v 

1 (u-1)H n (u) uH n _i(u) 

u ~ v (v-1 )H n (v) vH n _ r(u) 


(5.32) 


This formula generalize to arbitrary values of the vertex weights the result of paper [25] 
where an equivalent expression was derived in the case A = 7t/2 and rj = 7t/6, i.e., when 
a = b = c (the so-called ice point). 

As a final comment we would like to stress that formula flEUD implies that doubly 
refined weighted enumerations of ASMs can be easily obtained from the corresponding 
singly refined ones. In particular, the explicit expressions for the refined ^-enumerations, 
known for x — 1, 2, 3, can be just plugged into this formula to obtain the corresponding 
doubly refined ones. A discussion of the related one-point boundary correlation functions 
in application to the refined 1-, 2-, and 3-enumerations of ASMs can be found in [27]. 
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